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Abstract 

We demonstrate a method for exact determination of the quadratic curve of minimal energy and minimal curvature 
variation through three non-colinear points in the plane, including methods to determine the tangent vector and curvature 
at any point along the curve and an exact expression for the arc length of the curve between the first and last points. We 
then extended this to a novel method of selecting tangent vectors for use in constructing Hermite splines to smoothly 
O interpolate ordered sets of control points. Our results are spline curves of lower energy than that of many popular spline 
"implementations in most cases, which a series of examples demonstrate. 
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1. Introduction 

Interpolating splines have long been used to provide 
smooth curves through discrete sets of points. The ap- 
plications of such curves include interpolating data sets, 
computer-aided geometric design, typography, path plan- 
ning for robotics and CNC machining, highway and rail- 
way design, and in the design of computer games [l], 0, 
0, 0, 0, H 0, [H . Interpolating splines pass through each 
of an ordered set of defined control points, as opposed to 
estimating splines, which provide smooth curves that pass 
near their control points. 

There are many forms of interpolating spline, some 
more apt to a specific application than others. Two widely 
.used measurements of the quality of an interpolating spline 
are the energy (the elastic energy of the spline if it were 
an isomorphic flexible beam), and the curvature variation 
(how much does the curvature change over the length of 
the spline). Some general surveys, for example |9|, define 
additional criteria to judge interpolating spline quality, or 
fairness, including 

• extensionality (does adding a new data point on the 
spline change the spline?), 

• roundness (does the spline generate a circle given 
points lying on a circle?), 

• monotone curvature (do the curvature extrema fall 
on the control points?), and 



• locality vs. order (increasing order decreases ability 
to control spline locally). 

Some examples of interpolating splines currently in use or 
under research include 



• The family of parametric cubic curves [10|, UJj, in- 
cluding Hermite curves (of which Catmull-Rom [l2| 
and Kochanek-Bartels [13[ splines are special cases) 
and Bezier curves 

• Pythagorean-hodograph curves [3] 

• Minimal-energy curves [TBI 

• Minimal- variation curves [l6[ 

• Euler's spiral (or Conru spiral or clothoid curve) [ItJ 

• Log-aesthetic curves [l8[ 

Each of these have different strengths with respect to the 
criteria listed above. However, the family of spline curve 
that is most commonly used in design work are the para- 
metric curves, with cubics being the most common. 

There are an infinite number of such curves that can 
interpolate any given data set. In the case of Hermite 
curves, the choice of tangent vectors at each point selects 
from this family of curves, and the quality of the result- 
ing interpolating curve is based strongly on the method 
of choosing these tangent vectors. Before we discuss our 
method for selecting these tangent vectors, we briefly ex- 
plore parametric quadratic curves. 
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2. Quadratic parametric curves 

By parametric quadratic curves, we mean parametric 
curves represented by functions of the form 



r(t) = ait 2 + a 2 t + a 3 . 



(i) 



The curvature k, elastic energy E and curvature variation 
V of a parametric curve are given by 



k = 



xy — yx 



(x 2 + y 2 ) 3 / 2 ' 
/ K(t) 2 dt, 



V 



J m 2 



dt. 



(2) 

(3) 
(4) 



where dots indicate derivatives with respect to the curve 
parameter. The first and second derivatives of (pp) are 

r(t) = 2ait + a 2 and r(t) = 2a x , (5) 

and using a^ = (xi,yi),i G {1, 2, 3}, we have 

u\ o ^22/1 - X1V2 , a , 
K,(t)=2 o~7o" • (6) 

[(2^ + x 2 )2 + (2^ + ?/2 )2] 3 / 2 
Substituting (j6]) into (|3]), 

£ = 4(^22/1 - ^i2/2) 2 y [(2a?it + x 2 ) 2 + (2^/1 1 + ^ 2 ) 2 ] ~ 3 d£ . 
Evaluating this integral over the whole real line yields 

4 |£i2/2 - x 2 ?/i| 3 ' 
or in a coordinate-free form, 



E 



3tt 



4 ai x a 2 



3' 



(7) 



Then, differentiating (|2]) and using the fact that third 
derivatives vanish for a quadratic curve, 



> (xy — yx)(xx + \ 
1 (i 2 + y 2 ) 5 / 2 



(8) 



Applying ©, 



/i" = 12 



(xiy 2 - x 2 yi) (2 (x\ + y 2 ) t + xix 2 + 2/12/2) 



(2xit + x 2 ) 2 + (2 yi t + y 2 y 



5/2 



and so the curvature variation is 
V = 144 (xiy 2 - x 2 y x f 



I 



[2 (x\ ^y\)t^x 1 x 2 +yiy 2 \ 
{2x x t + x 2 f + (2 Vl t + y 2 f 



dt. 



Evaluating this over the real line gives the result 
45n(xl + yf) 2 



V 



16|#i2/2 - x 2 yi 



yf[(xi +yl) 2 \x\y2 -x 2 yi\ 2 

+ 7x 2 1 (x 1 -l) 2 (x 2 1 +y 1 y 2 ) 2 ] . 



However, we are more interested in minimizing V rather 
than obtaining its exact value, and so we introduce the 
following useful Lemma, 

Lemma 1. A quadratic curve of least energy also has least 
curvature variation. 

Proof. Consider a quadratic curve defined parametrically 
by y = at 2 + bt + c, with tangent vector given by y = 2at + b. 
The curvature at a point t is given by 

n = 2a [l + (2a£ + 6) 2 ]" 3/2 , 
with rate of change given by 

k = -12a 2 (2at + 6) [l + (2at + b) 2 ] " 5/2 . 

The total energy E and curvature variation V of the curve 
are given by 



E = Aa 2 / [1 + (2at + b) 2 ] 3 



(9) 



(2at + b) 2 [1 + (2at + b) 2 ] " 5 dt . (10) 

-00 



Integrating (0|) and HQ) , we obtain 



and V = — — lal 
16 1 1 



Therefore, the quadratic that minimizes \a\ will minimize 
both energy and curvature variation. 

3. Quadratic curve passing through three points 

Suppose we have three non-colinear points Pi, P2, P3 G 
R 2 . As pointed out in [l9|, these points can be interpolated 
by a parabola with r(0) = pi, r(l) = P3, and r(0.5) = P2 
(what Lachance and Schwartz call the "Cinci Parabola"). 
However, there is no reason to require that r(t) pass P2 
at t = 0.5. Rather, we wish to find the quadratic curve 
r(t) of least energy such that r(0) = pi, r(l) = P3, and 
r(T) = p 2 where < T < 1. Then, 

pi =r(0) = a 3 , 

p 2 =r(T) = aiT 2 + a 2 T + a 3 , and 
p 3 =r(l) = ai + a 2 + a 3 . 
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Solving for ai, a 2 , and a3, 

P2 - Pi - (P3 - Pl)T 



ai 

a 2 =P3 - Pi 

a 3 =Pi • 



rp2 _ 

P2 ~ Pi ~ (P3 ~ Pl)T 

rp2 _ J 1 



and 



(ii) 



We transform coordinates so pi lies at the origin, = 
Pi — Pi, i = 1,2,3, and then scale so |p 3 — pi| = 1, using 



P* 



P* 



IPs - Pi | 



Pi ~ Pi 

|P3 - Pi I 



1,2,3. 



We now use P3 to generate rotation matrix A that will 
rotate that point to (1,0), 



A 



PS X PSy 
-PSy P3 X 



Applying A to each point generates working points = 
Api. In this new coordinate frame, qi = 0, q3 = (1,0), 
and q2 y ^ 0. Equations ([TT]) are equally valid in this new 
frame, in which case they simplify somewhat, 



ai 



a 2 =q3 
a 3 =0. 



q2 - q3^ 
q2 - ^ 3 T 



2^2 



and 



(12) 



3.1. Minimum energy quadratic 
The energy, in terms of T, is 



E 



37ra |q 2 - q 3 T\ 4 
4|q 2 xq 3 | 3 T*-T 



To find the value of T for which the curve has minimal 
energy, we set, 



d [(q 2 - q 3 T) • (q 2 - q 3 T)] 2 



dT 



rp2 



0, 



Solving for T, recalling that |q3 1 
equation in T, 



1, generates a cubic 



(<?2 



|q 2 | 2 )T- 



102 l 



0. 



(13) 



When T = 0, the left-hand side of (fT3|) is positive definite, 
and when T = 1, the left-hand side reduces to 

~\ [(1 " <?2 J 2 + <Z2 y 2 ] , 

which is negative definite since q2 y ^ 0. Therefore, by 
continuity of ([T3]h we can be assured of having a root in 
the range < T < 1. Moreover, since ([T3]) goes to oo as 
T — >• oo and goes to — oo as T — oo, (fT3|) has three real 
roots, and we seek the middle root of the three. 



The roots of the cubic can be found using the cubic 
formula, 



T 2 



1 , M + C 
"2 + 2 ' 

i _ n + c - WHp- - 

"2 4 



, and (14) 



1 /u + C + *V3(A* - 



where 



[4(g 2:c -|q 2 | 2 )-3] 3 



/3 = 1 - 2(7 2;E , and 7 = 



27 



For example, suppose q 2 = (|, l) , in which case we would 
expect T = | by symmetry. In this case, /3 = and 
7 = —8, giving \i = iy^ and £ = —iy/2, so the roots are 



1 

2' 



T 2 



v^6 
2 ' 



and T 3 



v^6 
2 



and indeed, we obtain three real roots, one of which lies 
in (0, 1), that root falling at T = | as expected. 

Once T (the root between and 1) has been identified, 
we use (fTTj) to recover the a^ of the desired minimum en- 
ergy curve, then (pQ) to generate the curve. Reversing the 
coordinate transformations is not necessary since the value 
T is independent of coordinate system. A summary of the 
algorithm presented here, in more convenient pseudocode 
form, is included in the Appendix. 



4. Tangent vector at point p2 

To construct a Hermite spline through a sequence of 
points in which (pi, p 2 , P3) is a subsequence, we construct 
tangent vector r(T) at p 2 to the minimal-energy quadratic 
found above. Having solved for T, we have 



r(T) = (2T - 1) 



P2 - Pi - (P3 ~Pl)T 



rp2 _ J 1 



P3 - Pi 



For example, if pi = (0,0), p 2 = (^,1), and pi = 
(1,0), then r(T) = (1,0). Some examples of minimum- 
energy quadratic functions and the corresponding tangent 
vector are shown in Figure [TJ 

5. Arc length between pi and P3 

It may also be useful to know the arc length I of the 
curve between points pi and P3. The arc length is found 
by integrating the norm of the tangent vector over the 
interval. 
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Figure 1: Examples of ordered three-point sets and the corresponding 
minimum-energy quadratic. The tangent vector (derivative of the 
curve with respect to its parameter) at the center point is also shown. 



I = [ |2tai + a 2 |dt 
Jo 

= £ 2 ^P2-Pl-(P3-P1)T 



rp2 _ 



P3 - Pi 



P2 - Pi - (P3 - Pl)T 



rp2 _ J 1 

We define S3 = P3 — pi, and S2 = P2 — Pi, so 



dt. 



I 



S3 



dt 



J" _ r J 1 2 



[(2t-l)(Ts 3 -s 2 ) + (T-T 2 )s 3 ] 



dt. 



Then we note that for < T < 1, T - T 2 > 0, so 

1 r 1 

l= T=T*J |2t(^3-s 2 )-(T 2 S3-s 2 )| dt. 

Defining 1*1 = TS3 — s 2 and r 2 = T 2 ss — s 2 , and 6 as the 
angle between 1*1 and r 2 , 



I 



2 -4t|ri||r 2 |cose+ |r 2 | 2 dt. 



Evaluating this integral, we obtain 
1 



I = 



4|n|(r-r2) 



r 2 \ 2 cos<9 + (2|ri| - |r 2 |cos(9) p+ 



r 2 | sin ^log 



2|ri| - |r 2 |cos(9 + p 
|r 2 |(l-cosfl) 



6. A cubic Hermite spline implementation 

A Hermite cubic curve between any two consecutive 
points in an ordered point set is based on the locations of 
the points and on a tangent vector assigned at each point. 
The choice of tangent vectors at each vertex strongly af- 
fects the resulting spline. We now use the minimum- energy 
quadratic derived above to generate these tangent vectors, 
and compare the results to other methods. Note that the 
use of a single tangent vector at a control point for both 
adjoining spline segments assures G 1 continuity of the re- 
sulting curve. 




Kochanek-Bartels 



This paper's method 



Figure 2: Comparison of methods of choosing tangent vectors for 
cubic interpolating splines. The cardinal spline shown uses a tension 
parameter of 0.5. The Kochanek-Bartels spline uses a tension and 
continuity of and bias of 0.5. 

Figure [2] shows splines that result for several different 
methods of choosing of tangent vectors. If control points 
{pi, p 2 , p n } occur at parameter values {£1, t 2 , t n } 
along the spline, the methods of selecting tangent vectors 
{v 2 , V3, v n _i} shown in the Figure are listed here. 



6.1. Catmull-Rom spline: 



Pi+l — Pi-1 

U+i — U-i 



6.2. Cardinal spline: 



v< = (1 - r) 



Pi+1 ~ Pi-1 

U+i — U-i 
where r is a tension parameter, 



where p = y4|rip — 4|ri||r 2 | cos# + |r 2 | 2 . 
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6.3. Kochabek-B artels spline: 



_ (l-r)(l + /?)(!+ 7) , , 

v i — ~ (Pi - Pi-l) 



(l-r)(l-/3)(l- 7 ) 



(Pi+l - Pi) 



where r is a tension parameter, f3 is a bias parameter, and 
7 is a continuity parameter. 

6.4- A spline based on minimum- energy quadratics: 



H+i — H-i 



(2T- 1) 



Pi-i ~ (Pi+i - Pi-i)T 



rp2 _ J 1 



Pi+l - Pi-l | 



Where T is the root from (Q3]) that lies in < T < 1. 

We compare the elastic energies of spline segments be- 
tween the center two points of various four-point sets in 
the plane as shown in Figure OH 
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Point Set 1 

(0,0), (1,3), (2, 1), (3, 2) 



Point Set 2 

(0,0), (0,3), (3, 3), (3, 0) 



Point Set 3 



Point Set 4 



(0,0), (1,0), (2, 0), (3, 3) (0,0), (1,0), (2, 1), (3, 3) 

Figure 3: The four-point sets used to compare elastic energies and 
curvature variations of various spline curves. 

The elastic energies and curvature variations of our 
spline and standard splines with several parameter choices, 



as computed using E = J^ 3 k 2 dt, and V = J t ^ k 2 dt, are 
shown in Table [TJ 

The corresponding spline curves generated using mini- 
mum-energy quadratics are shown in Figure HI An exam- 
ple of a complete point set with the spline generated by 
the method outlined above is shown in Figure [5] 
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Set 1 



E: 


6.83 


13.46 


16.90 


71.53 


15.73 


9.41 


V: 


465 


1742 


2684 


41012 


2689 


943 


Set 2 














E: 


0.50 


0.50 


0.66 


3.65 


0.85 


0.85 


V: 


25.6 


25.6 


55 


2194 


15.1 


15.1 


Set 3 














E: 


0.49 


4.00 


3.69 


4.93 


1.57 


6.62 


V: 


14.6 


81.6 


78.6 


404 


16.4 


212 


Set 4 














E: 


0.45 


0.17 


0.19 


0.90 


0.69 


0.12 


V: 


3.0 


0.3 


1.8 


351 


10.2 


1.1 



Table 1: Energies and curvature variations of various splines, over 
the interval t<z <t <t$. 




Figure 4: The splines generated using the minimum-energy quadratic 
for the four test point sets. The energies in Table ^ are computed 
over the portion shown in red. 



7. Conclusion 

The method we have demonstrated in this paper to 
generate Hermite splines demonstrates significantly lower 
energies and total curvature variations than other com- 
monly used implementations for many (but not all) point 
sets. The calculations involved are more complex, to be 
sure, but these calculations need only be done during des- 
ign-time; calculating with or rendering the resulting spline 
curves can be accomplished using existing highly opti- 
mized Hermite curve routines. In interactive applications, 
this method can also be used to provide a reasonable de- 
fault set of tangent vectors at intermediate control points 
which act as a starting point for subsequent user manipu- 
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Figure 5: A randomly- generated ordered point set and the interpo- 
lating spline generated by building Hermite splines between succes- 
sive points where the tangent vectors at each point are those of the 
least-energy quadratic through the point and its two neighbors. 

lation. 

We acknowledge that there are point sets for which the 
splines generated by this method have higher energy and 
curvature variation than those generated by other methods 
such as Catmull-Rom. Having this method available, how- 
ever, would allow a designer to compute the energies for 
each potential spline implementation and select the opti- 
mal spline based on energy and curvature variation, which 
in many cases would be the spline whose tangent vectors 
are found using our method. 

Appendix A. Implementation algorithm 

This appendix gives an algorithm for Hermite spline 
generation based on the minimum-energy quadratic. 

Inputs: 

Four points p^ = (a^,^), i = 1,2,3,4, and a value of 
interpolation parameter t in the range [0,1], where t = 
corresponds to point P2 and t = 1 corresponds to point 

P3- 

Output: 

The interpolated point r(t) = (xt,yt) on the segment 
of the curve between p2 and P3. 

Algorithm: 

Procedure 1: To be performed when the set of control 
points changes. 



FOR EACH consecutive sequence of four points pi,p 2 ,P3,P4 

SET v 2 = C0MPUTE_TAN(pi,p 2 ,p 3 ) 

SET v 3 = C0MPUTE_TAN(p 2 ,P3,P4) 
END FOR 

FUNCTION C0MPUTE_TAN(q, r, s) 

// Shift first point to origin, scale so |s| = 1 

f = (r — a) /|s — a|, s = (s — a) /|s — a| 



// Rotate such that point s goes to (1,0) 

q = 0, r = (r x s x + r y s y , — r x s y + f y s x ) , s = (1, 0) 

// Compute the three roots of the cubic 

P = l-2r X9 7 =(4(f x -|f| 2 )-3) 3 /27, 1= ^=7 

01 = tan- 1 (V-T -/? 2 //?) /3 
\± r = I cos 4>i , \±i = I sin (f)\ 
fa = tan- 1 (-V-7-/3 2 //?) /3 
£ r = / cos 4>2 , Ci = I sin cf)2 

Ti = 1/2 + fi r + Cr/2 

T 2 = 1/2 - l/4[/x r + Cr + y/3(m - Ci)] 
T 3 = 1/2 - l/4[/x r + Cr - y/Sfai ~ Ci)} 

// Find the root between and 1 
IF Ti > AND Ti < 1 THEN 

T = Ti 

ELSEIF T 2 > AND T 2 < 1 THEN 

T = T 2 

ELSEIF T 3 > AND T 3 < 1 THEN 

T = T 3 
END IF 

// Construct the ai and a 2 coefficients 

ai = [(r-q)-r( S -q)]/(r 2 -T) 
a 2 = s — q — ai 

// Build and return the tangent vector v 
RETURN v = 2aiT + a 2 

Procedure 2: Hermite cubic spline interpolation, performed for each 
interpolation, with a value of t between and 1. 



r(t) = (2t 3 -3t 2 + l)p 2 + (-2t 3 + 3t 2 )p3 + (t 3 -2t 2 + t)v 2 + (t 3 -t 2 )v 3 
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